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ABSTRACT 

We examine the gravitational redshift of radiation emitted from within the potential of a cluster. 
Spectral lines from the intracluster medium (ICM) are redshifted in proportion to the emission- weighted 
mean potential along the line of sight, amounting to w 50 km/s at a radius of 100 kpc/h, for a cluster 
dispersion of 1200 km/s. We show that the relative redshifts of different ionization states of metals 
in the ICM provide a unique probe of the three-dimensional matter distribution. An examination of 
the reported peculiar velocities of cD galaxies in well studied Abell clusters reveals they are typically 
redshifted by an average of ~ +200 km/s. This can be achieved by gravity with the addition of a steep 
central potential associated with the cD galaxy. Note that in general gravitational redshifts cause a small 
overestimate of the recessional velocities of clusters by an average of ~ 20 km/s. 

Subject headings: galaxies: clusters: general — cosmology: observations — X-rays: galaxies — 
techniques: spectroscopic 



1. INTRODUCTION 

Recent N-body work has revealed that the central mass 
profile of cluster-sized dark matter halos in pure Cold Dark 
Matter (CDM) simulations are relatively flat inside a ra- 
dius of of ~ 400 kpc, asymptoting to r _1 (Navarro, Frenk, 
& White 1997, hereafter NFW; Moore et al. 1998). This 
prediction, although distinct, is not easily examined with 
lensing observations. The distortion of images by lensing is 
a powerful way of locating dark matter but is insensitive to 
the gradient of the distribution, particularly in the central 
region which is subject to a degeneracy (Schneider & Seitz 
1995; Kaiser 1995). Thus the flattening of the slope to r _1 
cannot be distinguished in practice from a singular isother- 
mal profile from lensing distortions alone (e.g. Clowe et al. 
1998). In principle, the magnifications of lensed sources 
are more sensitive to the projected mass profile, but their 
use requires knowledge of both the source redshifts and 
intrinsic unlensed properties to be reliable (Broadhurst, 
Taylor, & Peacock 1995; Taylor et al. 1998). Of course all 
lensing signals relate only to the projected mass distribu- 
tion and thus provide only indirect information along the 
line of sight. 

Here we evaluate a means of exploring the distribution 
of matter in clusters that complements and extends the 
information available from lensing. The gravitational po- 
tential of a cluster not only deflects light passing nearby, 
but imprints a gravitational redshift on light emitted by 
the cluster itself. The many X-ray emission lines in the 
metal-rich intracluster medium provide a suitable tracer 
of this redshift. Since the line emission is generated by 
a two body process, the velocity information is weighted 
by the square of the gas density along a given sight-line, 
and is quite sensitive to the mass profile in the central re- 
gion of the cluster. Furthermore, temperature gradients 
in the gas cause different lines to be emitted at different 
distances along a line of sight through the cluster, and 
may thus provide a three-dimensional probe of the matter 
profile. 



While X-ray spectroscopy has been difficult in the past, 
it is advancing to the point now where this effect may be 
marginally detected. In the near future, an order of mag- 
nitude increase in resolution and effective collecting area 
will be achieved with improved techniques. 

The structure of this paper is as follows. In §2 we calcu- 
late the form of the gravitational redshift for simple mass 
profiles. In §3 we discuss the properties of more realistic 
cluster profiles and the possibility of using multiple metal 
lines to study the distribution of matter along the line of 
sight. Finally in §4 we make some comments regarding 
existing optical data and prospects for the near future. 

2. GRAVITATIONAL REDSHIFT 

2.1. Simple Analytical Model 

The simplest useful model of a cluster is an isothermal 
sphere in which the density n(r) is proportional to r~ 2 for 
radii less than a truncation radius Rq, and n(r) — for 
r > Rq, so that the gravitational potential is given by 



( ) _ f 2cT 2 [ln(r/i?o) 



r < R 
r > Rq, 



(1) 



where a is the velocity dispersion. The gravitational red- 
shift of light from point sources such as galaxies located 
within the potential is given simply by 



Ai/(r) 



<f>(r) 



(2) 



where vq is the rest-frame emission frequency. Spec- 
tral line emission from the gas, however, is proportional 
to n e {r)n- lon (r) and hence assuming the metals are dis- 
tributed like the electrons we may simply weight the po- 
tential by n 2 (r) along lines of sight to compute the average 
redshift. In this case, the frequency of the emission line 
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will be shifted by an amount 

c 2 / fl(x) tfln((s2+I 2 )V2)2 

(3) 

where i is the distance along the line of sight and R(x) = 
For a singular isothermal sphere then, 
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where 6^ = arccos(x/i?o), L(6 X ) = — f Q * d8 ln(cos 9) is 
Lobachcvskiy's function, and the approximation is valid 
for small impact parameters. The velocity dispersion 
is only a weak function of cosmology (Eke, Cole, & 

1 /3 

Frenk 1996) and can be estimated as a — 730 Mi,f(l + 
z c ) 1 / 2 km/s, where Mu is the mass of the collapsed halo 
in units of 10 14 /i~ 4 solar masses and z c is the redshift of 
collapse. Finally, the relation between the total mass and 
the velocity dispersion (^jp = 2er 2 ) gives us the truncation 

radius R = 400 M 1 1 4 /3 (l + z c )- x kpc/h. 

In the top panel of Figure El we plot the redshift as a 
function of radius for three clusters with M14 = 1.5, 3, 
and 6, and z c = 0.5. Here we see that while the gravita- 
tional redshift of emission lines near the edge of the cluster 
is quite small, the r -4 weighting boosts the effect signif- 
icantly along lines of sight near the center of the cluster, 
reaching ~ 10% of the velocity dispersion within the inner 
10 to 50 kpc. 

2.2. Simple Numerical Model 

To better treat the interesting core region of the clus- 
ter, let us consider a simple nonsingular model (Binney 
& Tremaine 1987) derived by numerically integrating the 
equation of hydrostatic support of an isothermal gas 



d_ 

dr 



_ 2 d\np 
dr 



(5) 



where the density and radius have been rescaled such that 
p = p/po and r = r/ro, where po is the central density of 

the cluster and rp is the King radius, ro ~ 



4irGpo 



Inte- 



grating outwards from the boundary conditions p(0) = 1 
and jf = determines the density and gravitational po- 
tential at all radii, and we truncate the cluster at 10ro. 

Again we can relate mass and velocity dispersion and 
solve for ro as a function of mass and redshift, which gives 

r a = 38M^4 3 (1 + z c ) _1 kpc/h. We then use <j)(r) and n{r) 
to numerically integrate Eq. || to obtain the mean redshift 
of lines as a function of impact parameter x. In the cen- 
ter panel of Figure |l| we plot our results for three clusters 
of mass M14 = 1.5, 3, 6 and z c — 0.5, as the results for a 
nonsingular sphere in which the potential associated with a 
central cD galaxy has been simply added to the numerical 
solution, modeled as a point mass with 2% of the cluster 
mass. For reference, we also plot 4>(r) for the singular, 



nonsingular, and nonsingular plus point mass models in 
the bottom panel. 

In this figure we see that the singular and nonsingular 
isothermal models have much the same z(x) profiles out- 
side the core radius, but the presence of a core creates a 
sharp turnover at small radii, allowing these two poten- 
tials to be easily distinguished. Note also that the sharp 
increase in <f>(r) at very small radii in the core plus point 
mass model leads to only a ~ 10 km/s increase in the ob- 
served line redshifts, due to averaging along the line of 
sight. 

3. NFW MODEL AND TEMPERATURE GRADIENTS 

The models adopted in §2 were intended only as illus- 
trative, and by no means represent the full physics of the 
cluster dark matter and gas. The evolution of the dark 
matter in clusters is complex, many showing clear evidence 
for recent mergers (Zabludoff & Zaritsky 1995; Bird 1995; 
Bird, Davis, & Beers 1995; Bliton et al. 1998; Burns et al. 
2000). For many quiescent clusters, however, the radial 
profile of simulated clusters in CDM models are given by 
a universal density profile as shown in NFW: 



p oc 



r/r s (l + r/r s ) 2 



(6) 



where r s = ^-kpc/ hM^ 3 (1 + z c )'\ C « 3.6Q + z c ), and 
the profile is truncated at r/r s = C. Figure g illustrates 
one such distribution with M14 = 3 and z c = 0.5. Here 
we see that the steepness of the density profile and the 
flattening of the potential conspire to give a gravitational 
redshift that closely traces </>(r) and is very different from 
the singular isothermal model. Unlike the simple models 
considered in §2, however, the redshift estimate provided 
by Eq. || does not tell the full story, as it does not account 
for the presence of a temperature gradient. 

In general the gravitational redshift will not only be a 
function of gas density along the line of sight, but of the 
temperature profile as well. The ionization potential of 
each emission line will be strongly peaked about its ion- 
ization temperature (e.g. Zombeck 1990), and hence an 
emitted line will be confined to within a hollow shell in 
the cluster. Thus, with measurements of different X-ray 
emission lines, one can hope to simultaneously solve for 
the potential and the temperature profile, gaining three- 
dimensional information. This must be taken into account 
in a NFW potential and will be especially important for 
clusters containing cooling flows. In this case the tempera- 
ture gradient may exceed two orders of magnitude, so that 
for example the excitation of iron may range from FeX to 
FeXXVI (Fabian 1994; Churazov et al. 1998). 

A simple illustration of the additional information af- 
forded by temperature variations is given in Figure |^. Here 
the fractional ionization between two ionization states 
and Xi is allowed to vary linearly from lro to 3ro, 
in the potential given in §2.2, with M14 = 3 and both 
with and without a point mass. We interpret this model 
as a cooling-flow configuration in which the higher ion- 
ization state cannot exist within the core region of the 
cluster where the temperature drops much lower than its 
ionization potential. Here we see that the x beyond which 
Xi is completely ionized to the higher state is marked by 
an abrupt cutoff in line emission. The x at which Xi + \ 
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completely recombines to the lower state is marked by a 
leveling off in the redshift as a function of x, accompanied 
by a continuing increase in redshift as a function of x in 
the lines from the lower ionization state. Along lines of 
sight such that tq < x < 3ro, which at £ = pass through 
the partially ionized region, we see that the lines are offset 
from each other, providing information as to the width of 
the transition region along the line of sight. This can be 
compared to the width of the partially ionized region in 
the x direction to determine the sphericity of the cluster. 

4. OBSERVATIONS 
4.1. Redshift of Central Star Light 

While X-ray observations are the cleanest approach, the 
redshift of spectral features from starlight emitted from 
within the potential may also be studied as first suggested 
by Cappi (1995). Carefully studied massive clusters with 
central cD galaxies can be found in Sharpies, Ellis, & Gray 
(1988) and the sample of Zabludoff et al. (1993) and are 
listed below in Table |], excluding clusters with disper- 
sions below 700 km/s for which a significant redshift is 
not expected (Eq. ||). Interestingly, all these peculiar ve- 
locities but one are redshifted. In this table, the error 
listed in AV^d is the error in the measurement of the cD 
galaxy alone (ct c d) while the error listed in the following 
column is the error of the systemic velocity of the cluster 
(cr sys ). The total error on the peculiar cD velocity is thus 
<7gr = (er 2 D + o-gyj 1 / 2 . In some cases the clusters have well 
known cooling flows which we have noted below as the cen- 
tral potential may have been deepened by the action of a 
substantial flow over time. Here we use the values from the 
ROSAT study of Peres et al. (1998), quoting results ob- 
tained using the Position Sensitive Proportional Counter 
(PSPC) when multiple values were given. 

Substructure has been sought thoroughly for the above 
clusters and can always be identified at some level (Hill 
& Oegerle 1993, 1998; Oegerle & Hill 1994; Bird 1994; 
Kriessler & Beers 1997) but is not enough to account for 
the systematic redshift. Note we have excluded the cluster 
A85 for which significant clustering may be responsible for 
a large peculiar redshift of +1100 km/s (Bird 1994). 

Weighting the sum by the variance, we arrive at an av- 
erage AV c d of 260 ± 58 km/s. While this analysis is by no 
means exhaustive and does not attempt any correction for 
substructure, we may feel confident that the sign of this 
effect is real, namely a redshift, and the rough order of 
magnitude is not unreasonable. From Eq. || we can relate 
this redshift to the cluster potential. The error weighted 
ratio of AV^d and the velocity dispersion is (60 ± 20)cr 2 /c, 
or 340 ± 110 km/s for a — 1300 km/s, which, while some- 
what higher than our predicted redshifts, is nevertheless 
reasonable and may suggest that the form of the central 
potential is steeper than the slow logarithmic radial vari- 
ation of a purely isothermal potential. 

4.2. Future Prospects and Observational Issues 

While the systematic redshift of cD galaxies examined 
above provides a tantalizing suggestion of cluster grav- 

1 http://space.mit.edu/HETG/ 

2 http: / / astro.estec.esa.nl/XMM/ user / uhb / xmm_uhb.html 
3 http: / /constellation. gsfc.nasa.gov/design/resolution. html 
4 http: / /astro. estec.esa.nl/SA-general/Projects/XEUS/ 



itational redshifts, the possibility of subtle substructure 
makes studying optical emission from cluster galaxies dif- 
ficult to apply accurately. At best such studies are limited 
by the finite number of bright cluster members for mea- 
suring accurate velocities. 

Thus the X-ray spectroscopic approach outlined in §2 
and §3 will provide the most straightforward method for 
studying gravitational redshifts. Note that this requires 
both high spectral resolution to determine the shift of the 
line centroids, and fair spatial resolution to obtain radial 
information. Little information can be learned from a 
purely spectral analysis, however, as the overall cluster- 
integrated line profile is indistinguishable from a thermal 
profile with a slightly broadened width. 

The exploration of this effect may be marginally possi- 
ble with the current generation of X-ray satellites. The 
High-Energy Transmission Grating on the Chandra Satel- 
lite, launched in July 1999, is already reaching a resolv- 
ing power of E/AE « 1,000 = 300 km/s in the 1 to 2 
keV range, 1 and the Reflection Grating Spectrometer on 
the XMM satellite launched in December 1999 provides 
E/AE values from 200 to 800 in the energy range 0.35 
to 2.5 keV for point sources, although the slitless design 
of this instrument reduces the spectral resolution of ex- 
tended sources such as clusters. 2 Typically, using mul- 
tiple emission lines, one can achieve resolutions of 10% of 
these values, suggesting that investigations of gravitational 
redshifts in rich clusters may be possible but challenging. 
Note that the thermal velocity of the ions is smaller than 
a by a factor of yj m p /m.i on . 

The next generation of satellites will significantly in- 
crease these capabilities. The Spectroscopy X-ray Tele- 
scope (SXT) on the proposed Constellation X team 
of X-ray telescopes plans to use both a quantum mi- 
crocalorimter with 2 eV resolution and a set of reflection 
gratings below 2 keV, to obtain velocity sensitives bet- 
ter than 100 km/s below 1 keV and of the order of 20 
km/s at higher energies. 3 The narrow field detectors on 
the proposed XEUS observatory in space promise simi- 
lar resolutions. 4 With such sensitivities, both core-profile 
studies such as shown in Figure [j] and multiple-line, three- 
dimensional studies as described in §3 will be well within 
reach. 

Some observational issues are difficult to predict. As line 
emission is proportional to n e (r)ni on (r) , metallicity gradi- 
ents can bias gravitational redshift measurements. Such 
gradients have already been discovered in the Centaurus 
Cluster (Fukazawa et al. 1994) and the cluster AWM 7 
(Ezawa et al. 1997), and will undoubtedly be a subject of 
investigation for the current generation of X-ray satellites. 
As in some cases the metallicity can vary by a factor of two 
over the cluster, these gradients must be well understood 
before our approach is applied. 

Large-scale bulk motion of gas should also be detected 
from recent group infall in some clusters. The disturbance 
of the intracluster medium by gas infall has been claimed 
to be significant over a period of 3 Gyrs in favorable condi- 
tions, before shocks are thermalised and bulk flow through 
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the core ceases (Roettiger, Burns, & Loken 1996) If infall is 
common then the relatively small perturbation to the line 
centroid by gravity proposed here may be made noisy by 
the continuous process of cluster interaction, or perhaps 
in some cases relegated to a correction applied to a larger 
dynamical disturbance. Finally the rclativistic tempera- 
tures in the cluster gas will cause a small shift in X-ray 
line emission by a factor of 7 — 1 kT(r) due to rel- 
ativistic beaming and the transverse Doppler effect, but 
fortunately this is smaller than the gravitational redshift 
by a factor of m p /m- lon . 



These and other issues will await the new and next gen- 
erations of X-ray telescopes. It is clear, however, that 
whatever the experimental challenges, a careful measure- 
ment of X-ray emission lines over the surface of galaxy 
clusters will be a useful probe of cluster structure and for- 
mation. 
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Fig. 1. — Top: Redshift as a function of impact parameter for a singular isothermal profile with M14 = 1.5, cr = 1000 (solid), M14 = 3, 
cr = 1300 (dotted), M14 = 6, cr = 1600 (dashed). Center: Redshift as a function of impact parameter for a nonsingular isothermal profile 
with masses and velocity dispersions as in the upper panel. In each pair of lines the lower line represents the redshift purely from the profile, 
while the upper line includes an additional contribution from a cD galaxy modeled as a point mass with 2% of the cluster mass. Bottom: 
Gravitational potential for three different density profiles. The solid lines represent the nonsingular isothermal sphere, the dotted lines show 
the singular isothermal model, and the dashed lines show a nonsingular model in which a cD galaxy has been added. In all cases M14 = 3, 
o- = 1300, and z c = 0.5 (see §2.1). 




Fig. 2. — NFW model for a M14 = 3, z c = 0.5 cluster. The upper panels show the density profile, cluster potential, and temperature 
(velocity dispersion) as a function of radius. The lower panel shows the "naive" gravitational redshift as a function of impact parameter as 
predicted by Eq. H. 
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Fig. 3. — Gravitational rcdshift for the Mi 4 = 3 model in Figure 2 in which the effects of temperature variations are modeled simply by 
taking the fractional ionization between two ionization states to vary linearly from lro to 3ro. The lower ionization state is given by the solid 
lines and the higher by the dotted lines. Top: Observed gravitational redshift of the lines from the two ions as a function of impact parameter 
of the line of sight. Bottom: The fractional number densities of each ionization state as a function of radial distance. 



Table 1 
Anomalous cD Redshifts 



Cluster 




Dispersion 


AV cD 


<^sys 


Cooling 


Reference 






(km/s) 


(km/s) 


(km/s) 


Flow 

(M /yr) 




A193A 


32 


760 ± 100 


240 ± 150 


150 




Zabludoff et al. (1993) 


A426 


114 


1277 ± 80 


-106 ± 124 


122 


536+24 


Zabludoff et al. (1993) 


A1644A 


82 


940 ± 70 


35 ± 150 


113 


H + 5° 


Zabludoff et al. (1993) 


A1651 


29 


965 ± 130 


210 ± 200 


195 


138t£ 


Zabludoff, Huchra, & Geller (1990) 


A1656A 


270 


1140 ± 50 


390 ± 80 


64 




Zabludoff et al. (1993) 


A 1795 


150 


773 ± 80 


360 ± 150 


124 


38lt£ 


Hill et al. (1988) 


A2199A 


68 


820 ± 70 


230 ± 110 


104 


154+f 


Zabludoff et al. (1993) 


A2690 


303 


880 ± 40 


440 ± 150 


84 




Sharpies, Ellis, & Gray (1988) 



